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Abstract. For a polynomial P, we consider the sequence of iterated integrals 
of lnP(a;). This sequence is expressed in terms of the zeros of P{x). In the 
special case of ln(l + z^), arithmetic properties of certain coefficients arising 
are described. Similar observations are made for ln(l + x'^). 



1. Introduction 

The evaluation of integrals, a subject that had an important role in the 19*'' 
century, has been given a new life with the development of symbolic mathematics 
software such as Mathematica or Maple. The question of indefinite integrals was 
provided with an algorithmic approach beginning with work of J. Liouville [8j dis- 
cussed in detail in Chapter IX of Lutzen [9]. A more modern treatment can be 
found in Ritt (21], R. H. Risch [HI [20], and M. Bronstein [3]. 

The absence of a complete algorithmic solution to the problem of evaluation of 
definite integrals justifies the validity of tables of integrals such as [H H] [TB] . These 
collections have not been superseded, yet, by the software mentioned above. 

The point of view illustrated in this paper is that the quest for evaluation of 
definite integrals may take the reader to unexpected parts of mathematics. This 
has been described by one of the authors in [Ml [15] . The goal here is to consider 
the sequence of iterated integrals of a function fo{x), defined by 

(1.1) fnix) = / fn-l{t)dt if n> 1. 

Jo 

This formula carries the implicit normalization /„(0) = for n > 1. 
A classical formula for the iterated integrals is given by 

(1-2) fnix) = - I hit) ix - tr-' dt. 

Expanding the kernel (x — i)"^^ gives /„ in terms of the moments 

(1.3) Mj{x) = / t^foit)dt 

as 







(1.4) fnix) = ^(-1)^—^- -M,ix). 

^ j!(n-l-j)! 
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The work presented here deals with the sequence starting at fo{x) ~ ln(l + x^). 
The main observation is that the closed-form expression of the iterated integrals 
contains a pure polynomial term and a linear combination of transcendental func- 
tions with polynomial coefficients. Some arithmetical properties of the pure poly- 
nomial term are described. 

2. The iterated integral of ln(l -t- x) 

The iterated integral of fo{x) — ln(l -I- x) was described in ^13j. This sequence 
has the form 

(2.1) /„(a;) = A„,i(x) +S„4(a;)ln(l-f a;) 
where 

(2.2) A„,(.) = -1 x: Q (H. - i/„-.) "'^"V^""' ' 

' fc=l ^ ^ ' k=l 

BnAx) - ^{i + xy\ 

where H„ ~l~'"i~^'''+n harmonic number. 

The expression for Bn,i{x) is easily guessed from the symbolic computation of 
the first few values. The corresponding closed form for A„.i(x) was more difficult to 
find experimentally. Its study began with the analysis of its denominators, denoted 
here by a„,i. The fact that the ratio 

(2.3) /3„,i := 
satisfies 

\ p if n is a power of the prime p 

I 1 otherwise 



(2.4) /?„ 



was the critical observation in obtaining the closed form An.i{x) given in ()2.2p . We 
recognize f3n,i as e^^^"), where 



(2.5) A(n) = 

is the von Mangoldt function. This yields 



Inp if n is a power of the prime p 
otherwise 



J=2 j=2 



A(n) ^ lcm(l,...,n) 
lcm(l, . . . , n — 1) 



and the relation 
(2.6) 
shows that 

(2.7) a„^i — nl lcm(l, . . . , n). 

Note 2.1. The harmonic number Hn appearing in (|2.2p has challenging arithmeti- 
cal properties. Written in reduced form as 

(2.8) 




the denominator D„ divides the least common muhiple i„ :— lcm(l, 2, . . . , n). The 
complexity of the ratio Ln/Dn can be seen in Figure [TJ It has been conjectured [3 
page 304] that £)„ = _L„ for infinitely many values of n. 

The expressions for An.i{x) and Bn^i{x) can also be derived from (11.41) . Letting 
fo{x) = ln(l + x) yields 

(2.9) /„(,)^^^(_i).____ tnnil + t)dt. 

Integration by parts gives 

r xi+^ln{l + x) 1 rp+^dt 

(2.10) / tnYY{l + t)dt= ^ ' ' 



J + l J + lJo 1 + i 

Replacing in (|2.9p shows that the contribution of the first term reduces simply to 
x" ln(l + x). Therefore 

(2.U) /„(., = 1." ,„(i + .) + 1 p-,y I' fii. 

It remains to provide a closed form for the integrals 

(2.12) '^■■-llT-t''- 
These can be produced by elementary methods by writing 

(2.13) j^^P-i-iy til. 
^ ' 1+t l+t 1+t 
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Replacing in (|2.11|) gives 



fn{x) = la;"ln(l + x) 



+ 



i!^Vj7 Jo 1 + 



- {-ly 
t + i 



dt 



The first and last line add up to (x + 1)" ln(l + x)/n\, which yields the closed-form 
expression for Bn^i{x). Expanding the quotient in the second line produces 



1 " 

(2-14) ^E(-l) 



E 

r=0 



The double sum can be written as 

71 n—j 



(2.15) 



-yy 

j=0 r=l 



,j+r _ 



n-1 



j=o 



11-] 



1 " 

-y 



E 



n v-ir 

a — r 



The expression for An^i{x) now follows from the identity 

n \ {-ly 



(2.16) 



E 



a — r J r 



[Hn — Hn-a] 



An equivalent form, with to = n — a, is given by 
(2.17) 



U{a) '■= E ~ — /L-4-i'''^ ~ Hm+a - Hjn- 



To establish this identity, we employ the WZ method [TB]. Define the pair of 
functions 



(2.18) 



Fir,a)^ ^ and G(r,a) 



(-l)'^C-i) 



(m + a + l)("\tV)' 



r ("+'■) 

It can be easily checked that 

(2.19) F(r, a + 1) - F(r, a) = G'(r + 1, a) - G(r, a). 
Summing both sides of this equation over r, from 1 to a + 1, leads to 

(2.20) U{a + l)-U{a) = 



1 



TO + a + 1 

Now sum this identity over a, from 1 to fc — 1, to obtain 

k — l m+k 



c/(fc) - L/(i) = E — ^ = E - 

■■^ TO + a + 1 „ r 



(2.21) 

a— 1 ' ' a— m+2 

Combining this with the initial condition U{\) ~ gives the result. 



m+l- 
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3. The method of roots 
The iterated integrals of the function /o(a;) = lnP(2;) for a general polynomial 

m 

(3.1) P{x) = l[{x + z,) 

i=i 

are now expressed in terms of the roots Zj using an explicit expression for the 
iterated integrals of fo{x) = ln{x + a). 

Theorem 3.1. The iterated integral of fo{x) = ln(x + a) is given by 

(3.2) U{x) = :V ^ , ^ ^ -. lna+^ ^ln(x + a). 

k=\ 

Proof. A symbolic calculation of the first few values suggests the ansatz /„(a;) = 
Sn{x) + Tn{x) hi{x + a) for some polynomials 5*^, T„. The relation = /n-i and 
the form of S'„,r„ given in (j3.2p give the result by induction. □ 

In the special case P{x) = 1 + x^ , the previous result can be made more explicit. 

Theorem 3.2. Let a ^ u + iv he a root of 1 + = 0. Then the contribution of 
a and a ~ u — iv to the iterated integral of ln(l + x^) is given by 

1 " „fc 



n! ^-^ k 

k=l 



-i- [{x + a)" — {x + a)"] arctan 



inl V 1 + ux 

ln[(l + uxf + v^x 



(x + a)" + (x + a)" , . ,^ , , 2 , 2 21 



2nl 

Proof. First observe that ln(a:; + a) — In a = ln(aa:; + 1); hence for fo{x) = ln{x + a) 
Theorem 13. II takes the form 

(3.3) fn{x) = V ^ — + — Ina + ^ ln(aa; + 1). 

n\ ^ — ' k n\ n\ 

k=l 

Since 

ln(ax + 1) = In \ax + 1| + z Arg(aa:; + 1), 
ln(aa; + 1) = In \ax + 1\ — i Arg(aa:; + 1), 

and Ina + Ina — 21n|a| = 0, it follows that the total contribution of a and a is 
given by 

1 V-^' , .^n-k , , I , [(x + a)" ln(ax + 1) + (x + a)" ln(ax + 1)] 



n! ^ k 

fc=i 



1 a;'' ,- , ^,,1. , _^„_^.^ [(a; + a)" + (a; + a)'M 
= -^2^ — [(x + a)" ^ + {x + aY ^]+- '—^ ^ln|ax + l| 

■ fc=i 

i [(x + a)" - (x + a)"] ^ , ^ , 
Lv )_ — \ Arg(ax + 1). 

n\ 

The stated result comes from expressing the logarithmic terms in their real and 
imaginary parts. □ 
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Corollary 3.3. Let n G N. Then 



r+k/,, \n-k 1 a; {x + a) 

Proof. Integrate both sides of the identity in Theorem 13.11 and use the relation 
fn-i — fn to obtain the result inductively. □ 

Note 3.4. The identity in Corollary 13.31 can be expressed in terms of the function 

" 1 

(3.4) $„(a;,a) i^^-x'^lx + ar-^- 
in the form 

(3.5) r $„(*, a) dt = „(.T, a) + , [x''+' + a"+i - {x + a)"+i] . 
Jo n + i (n + ij 

The function $„(a;, a) admits the hypergeometric representation 

^ r ^ 2;"+^ ^/1,1 + n X \ . f a \ 

$„ x,a^- : ;— — - x + a "In — — . 

(n + l)(a; + a) \/ + n x + a/ \x + a J 



With this representation, the identity in Corollarv l3.3l now becomes 



i \ ^ / 1, 1 + n ,\ at 1 I X ^ 



1-tl ^ 'V 2 + n '7l-i n+lVl-x 



/ 1, 1 + n- 1 1 
2^^i( 2 + n 



4. The iterated integral of ln(l + x'^) 
In this section we consider the iterated integral of fo{x) = ln(l + x'^) defined by 



(4.1) fnix) = / fn-l{t)dt. 

Jo 

The first few examples, given by 

fi{x) = — 2.T + 2 arctana; + x ln(l + x'^) 

/^{x) = + 2x arctanx + ^{x'^ - 1) ln(l + x'^) 

fsix) = + 1^ + (^^ ^ i) arctanx + [^x^ - ^x) ln(l + x^) 

suggest the form 

(4.2) fn{x) = An^2{x) + Bn,2{x) arctanx + Cn,2{x) ln(l + x^) 



for some polynomials v4„_2, Sn,2, C'„.2. Theorem 13.21 can be employed to obtain a 
closed form for these polynomials. It follows that fn{x) satisfies 

n J. 

(4.3) n\Ux) = - ^ ^ [(.X + ^r-^ + (x - z)"'^'] 
fe=i 

- i [(x + i)" - (a; - «)"] arctanx + ^ [(^^ + «)" + (a; - i)"] ln(l + x^). 
The expressions for A„_2, 5n,2, C'n,2 may be read from here. 
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4.1. Recurrences. The polynomials An^2, Bn.2, Cn,2 can also be found as solutions 
to certain recurrences. Differentiation of ()4.1|) yields fn{x) — fn-i{x). It is easy to 
check that this relation, with the initial conditions /„(0) = and /o(a;) = ln(l + a;^), 
is equivalent to (|4.1I) . Replacing the ansatz (|4.2p produces 



A'nA^) + B'^^x) arctanx + ^^^^ + C'^^x) ln(l + x') + ^,2(2;)^^ 

1 + ' 1 + a;^ 

= ^1-1,2(2^) + Bn-i^2ix) arctanx + C„_i^2(a;) ln(l + x'^). 
A natural linear independence assumption yields the system of recurrences 



(4.4) Bl^^x 
Bq^x 

(4.5) CUx 

Co,2ix 

(4.6) A'^^x 
AoM.x 



Bn^l^x) 



C„-i,2(a;) 
1 

- / X Bn,2{.x) + 2xCn,2{x) 
^-1,2(2:) 

0. 



Note 4.1. The definition (|4.ip determines completely the function fn{x). In par- 
ticular, given the form ()4.2p . the polynomials A„^2,Bn,2 and Cn,2 are uniquely 
specified. Observe however that the recurrence (|4.4p does not determine Bn^x) 
uniquely. At each step, there is a constant of integration to be determined. In order 
to address this ambiguity, the first few values of -6^.2(0) are determined empirically^ 
and the condition 



(4.7) B„,2(0) 



( 11—1 
^l2(-l)~/n! if n is odd 

if n is even 



is added to the recurrence (|4.4p . The polynomials B^^x) are now uniquely deter- 
mined. Similarly, the initial condition 



(4.8) C„,2(0) = 



(—1) 2 ln\ if n is even 
if n is odd 



adjoined to (|4.5p . determines C„ 2- The initial condition imposed on A„ 2 is simply 

A.,2(0) = 0. ' 

The recurrence (|4.4p is then employed to produce a list of the first few values of 
Bn,2{x). These are then used to guess the closed- form expression for this family. 
The same is true for Cn,2{x). 

Proposition 4.2. The recurrence (|4.4p and the (heuristic) initial condition (|4.7p 
yield 

n-l 

3=0 ^ ■' 



= -l[(x + z)"-(x-z)^ 
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Similarly, the polynomial Cn,2 is given by 

L"J 

(4.10) ^».2W-^E(-1V(2")-""'^ 

2 n\ 

In particular, the degree of Bn.2 is n — 1, and the degree of Cn,2 is n. 

Proof. This follows directly from the recurrences (|4.4p and (|4.5p . □ 

Corollary 4.3. The recurrence for An.2 can be written as 

(4.11) A;_2(.t) = A„_i,2(:e) --.[{x + z)"-i + (x - z)"-i] . 
In particular, the degree of An ^2 is n. 

Proof. Simply replace the explicit expressions for Bn,2 and Cn,2 in the recur- 
rence □ 

4.2. Trigonometric forms. A trigonometric form of the polynomials Bn,2 and 
Cn.2 is established next. 

Proposition 4.4. The polynomials i3„.2 and C„^2 are given by 

2 

Bn2(x) = — (x^ + 1)"/^ sin(narccotx) 



Cn2{x) = ^(x^ + 1)"/^ cos(narccotx). 



In particular, 

(4.12) ^2^44 = - cot(n arccot x). 

Bn.2(.x) 2 

Proof. The polar form 



(4.13) a; + i = vaS^ + l [cos(arccot x) + i sin(arccot x)] 
produces 

(4.14) {x + iy^ = (x^ + 1)"/^ [cos(n arccot x) + i sm(n arccot x)] . 

A similar expression for (x -~ i)" gives the result. □ 
Proof. A second proof follows from the Taylor series 



sin(zarctan^) ^ >^ (-l)''(^)2fc+i 2fc+i 
^ ' (l + t2)-/2 A. (2fc + l)' 

and 



cos(z arctani) (— l)'^(z) 



2k 2k 



(l+t2)-/2 "A. (2/fc)! 

where (z)„ denotes the Pochhammer symbol. These series were established in [5] 
in the context of integrals related to the Hurwitz zeta function. 
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Indeed, the formula for Bn,2{x) comes from ()4.15p replacing t hy 1/x and z by 
— n to obtain 

(4.17) sin(narccotx) (x^ + 1)"/^ = -a;" i^il^Mz^x-^'^-i. 



The result (|4.9p now follows from the identity 
(4.18) (-n)2fc+i 



-n\/{n-2k-iy. if2fc+l<n 
otherwise. 



A similar argument gives the form of Cn.2{x) in (I4.10p . □ 
Note 4.5. The rational function i?„ that gives 
(4.19) cot(n6') = i?„(cot6l) 

appears in (|4.12p in the form 

4.20 Rn{x) = ; ; . 

This rational function plays a crucial role in the development of rational Landen 
transformations [lOj . These are transformations of the coefficients of a rational 
integrand that preserve the value of a definite integral. For example, the map 

a ^ a ((a + icf - 35^) /A 

h ^ b{3{a-cf -b^) /A 

c ^-> c((3a + c)2 -362)/A, 

where A = (3a + c){a + 3c) — b^, preserves the value of 

dx 27r 



ax^ + ox + c v4ac — b^ 

The reader will find in 12 a survey of this type of transformation and [11) the 
example given above. The reason for the appearance of i?„(x) in the current context 
remains to be clarified. 

4.3. An automatic derivation of a recurrence for A„_2- The formula (|1.2p for 
for the iterated integral can be used in the context of computer algebra methods. 
In the case discussed here, the integral 

(4.22) In{x)^ . / {x-tr-Hn{l + t^)dt 

[n-iy- Jo 

gives the desired iterated integrals of ln(l + x^) for n> 1. 

A standard application of the holonomic systems approach, as implemented 
in the Mathematica package HolonomicFunctions [6|, yields a recurrence in n 
for (|4.22p . The reader will find in 7 a description of the use of this package in the 
evaluation of definite integrals. The recurrence 

(4.23) n2(n-l)/„(x) 

= x(3n - 2){n - 1) /„_i(x) - (Srix^ - Ax"^ + n) In-2{x) 

+ x(x2 + l)/„_3(x) 
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is delivered immediately by the package. Using the linear independence of arctana; 
and ln(l + x^), it follows that each of the sequences An,2, -Bn.2, and Cn,2 must 
also satisfy the recurrence (j4.23p . Symbolic methods for solving recurrences are 
employed next to produce the explicit expressions for An^2, 5n,2, and Cn,2 given 
above. 

Petkovsek's algorithm Hyper '17| (as implemented in the Mathematica package 
Hyper, for example) computes a basis of hypergeometric solutions of a linear re- 
currence with polynomial coefficients. Given (|4.23l) as input, it outputs the two 
solutions (x + and {x — i)'^/n\. The initial values are used to obtain the 

correct linear combinations of these solutions. This produces the expressions for 
Bn,2{x) and Cn,2{x) given in Proposition 1121 

However, the third solution is not hypergeometric and it will give the polynomi- 
als An,2{x). It can be found by Schneider's Mathematica package Sigma 22 : 

AnA^) ^-Ax {{x + ^r -{x~ iT) + -\ , 

with the initial values 

^o,2(a;)=0, Ai,2(a;) = -2a;, ^2,2(0;) = -|x^ 
In summary: 

Theorem 4.6. Define a/, = k{k — 1) for /c > 2 and ai = —1. The polynomial 
An,2{x) introduced in (14. 2p is given by 

(4.24) A„.2(x) = V — [(a; + 1)"-'^+^ ~ix- *)"-^+i] . 

k=\ 

This can be written as 

(4.25) A„,2(x) = ^ V ^"~^^^^' x^-i3„-fc+i,2(x). 

Note that the expression for A„^2 given before is equivalent to the forms appear- 
ing in Theorem 14.61 

Note 4.7. Similar procedures applied to the case of ln(l -)- x) yield the evaluation 
given in (j^ . 

5. Arithmetical properties 

In this section we discuss arithmetical properties of the polynomials -B„,2 and 
An,2- The explicit formula for i?„^2 produces some elementary results. 

Proposition 5.1. Let m,n e N such that m divides n. Then i?m,2(a;) divides 
Bn,2ix) as polynomials in Q[x]. 

Proof. This follows directly from (|4.9p and the divisibility of a" — 6" by a™ — . □ 

For odd n, the quotient of -B2n,2(a;) by B„^2ix) admits a simple expression. 
Proposition 5.2. Let n E N. Define 
(5.1) B:^,{x)^x''-^''"'-B,,^2il/x). 
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2n\ n-1 

B2n.2{x) = (-1) 2 xBnMx)B^^2(x)- 



Then, for n odd, 

In particular, the sequence of coefficients in B2n{x) is pahndromic. 
Proof. The proof is elementary. Observe that 



in\ 
1 

ixnl 

■n-1 



[(1 



I 

X 

(I - ixY 



[(x-^r-i-lnx + ^r] 



It follows that 



(5.3) 



K,2i^) 



n-1 

{-!) — 



[{x + 0" + ix~ lY' 



and the result now follows directly. □ 

The explicit expression (j4.24l) for the polynomial An^2 can be written in terms 
of the polynomials 

(5.4) = + ir - (x - tr 

as 



(5.5) 



An,2{x) 



Xipn{x) 



k=2 



k{k-l) 



The polynomial An^2 is of degree n and has rational coefficients. 

By analogy with the properties of denominators of v4„ i {x) mentioned in Section[2] 
and discussed at greater length in |13| . we now study the denominators ^,1^2(2;) from 
an arithmetic point of view. The first result is elementary. 

Proposition 5.3. Let 

(5.6) an,2 ■— denominator of An^2{x). 

Then a„.2 divides n\ lcm(l, 2, . . . , n). 

Proof. The result follows from (15. 5p and the fact that the polynomials (pm{x) have 
integer coefficients. □ 

As in (|2.3p . it is useful to consider the ratio 

On, 2 



(5.7) 



Pn,2 ■■ = 



nan-1,2 

Symbolic computations suggest the following. 

Conjecture 5.4. The sequence /3„,2 is given by 

p if 71 = p'' for some prime p and r E N and n ^ 2 ■ 3^' 

i if n = 2 • 3" for some m £ N 

(5.8) /3n,2 = S 3p if n = 2 • 3" + 1 and n = p'' for some m,r eN 

3 if n = 2 • 3'" + 1 for some m G N and n 7^ 

1 otherwise. 
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The formulation of this conjecture directly in terms of the denominators of 
An,2{x) is as follows. 

Conjecture 5.5. The denominator an,2 of An^2{x) is given by 

{1 if n = 1 

nl lcm(l, 2, . . . , n)/6 if n = 2 • 3™ for some m > 1 
n! lcm(l, 2, . . . , n)/2 otherwise. 

This conjecture shows that the cancellations produced by the polynomials ipm (x) 
in (15.51) have an arithmetical nature. 



Proof that Conjecture \5.5\ implies Conjecture \5.4\ Assume that (15.91) holds for n > 
1. If n = 2 • 3™, then a„_2 contains one fewer power of 3 than nan~-i,2- li n = 
2 • 3™ + 1, then a„ 2 contains one more power of 3 than nan-.i.2- li n = is a 
prime power, then q;„_2 contains one more power of p than nan-i.2- Otherwise 
each prime appears the same number of times in a„.2 a-nd T^ctn- 1,2- D 

The first reduction is obtained by expanding the inner sum in (j5.5p . Define 



L"/2J-1 

(5.10) Gn{x)=-2i J2 (-1)* 

k=Q 

Proposition 5.6. We have 



n-l ^ / 

y ^ f ' 

jjk+i in-j){n-j + l)\2k + l 



^n-2k 



(5.11) A^,2{x) = -\x {{x + ir ~{x~ iD + G,,{x)] . 

Proof. Expanding the terms {x + i)"-'=+i and (x — i)"-'=+i in the expression for 
An,2{x) yields the sum 

1 n+l-j i / -x 

^ (n-j)(?i-j + l)^ 

so only odd k contribute to it. Reversing the order of summation gives the result. 

□ 

The next result compares the denominator a„^2 of An ^2 (x) and the denominator 
of G„, denoted by 7„. 

Corollary 5.7. For n > 1, the denominators an,2 and jn satisfy 
(5.13) an,2 = n!7„. 

We now rephrase Conjecture 15.51 as the following. 
Conjecture 5.8. For n > 2, 

{lcm(l, 2, . . . , n)/6 if n 2 • 3™ for some m > 1 
lcm(l, 2, . . . , n)/2 otherwise. 

The next theorem establishes part of this conjecture, namely the exceptional role 
that the prime p = 3 plays. The proof employs the notation 

(5-15) 9n.k{j) = 



(5.14) 7„ 



in-j){n-j + l) V2fc + 1 
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SO that 



L"/2J-1 



(5.16) Gn{x)^-2i (-l)^-Vfea;"-2fc 



n~l n— 1 — 2/c -, / /) 

1 I n — t 



with 

(5.17) V.:= ^ 5„,.(j-)= E ^(^+l)V2fc + l 
Therefore, for n > 2, 

(5.18) 7n = 2 ■ Icm {denominator of /i„^fe : < A: < [n/2j — 1} . 

Let Vp{n) be the exponent of the highest power of p dividing n — the p-adic 
valuation ofn. The denominators in the terms forming the sum hn,k are consecutive 
integers bounded by n. Therefore 

(5.19) i^sijn) < z^3(lcm(l,2,...,n)). 
In fact we can estabUsh vai^n) precisely. 
Theorem 5.9. The 3-adic valuation of 7„ is given by 

z^3(lcm(l,2, . . . ,n)) — 1 if n = 2 • 3™ for some m > 1 
i^3(lcm(l, 2, . . . , n)) otherwise. 

Proof. The analysis is divided into two cases. 

Case 1. Assume that n = 2 • 3™. We show that lysijn) — m — 1. 

The bound (|5.19p shows that 2^3(7^) < m. 
Claim: I's^'^n) 7^ tti- To prove this, the coefficient 

is written as 

(5.21) hn,k = Siin,k) + S2in,k) 

where Si{n, k) is the sum of all the terms in /i„ /j with a denominator divisible by 
3™ and S2{n,k) contains the remaining terms. This is the highest possible power 
of 3 that appears in the denominator of hn^k- 

It is now shown that the denominator of the sum Si(n,k) is never divisible by 
3™. 

Step 1. The sum Si{n,k) contains at most two terms. 

Proof. The index £ satisfies £ < 2 • 3" - 1 - 2fc < 2 • 3". The only choices of £ that 
produce denominators divisible by 3™ are i — 3™, 3™ — 1 and £ = 2 • 3™ — 1. The 
term corresponding to this last choice is (2-3'" -\) ■2-3'" (2fc+i) ' only occurs for 
k — 0. In this situation, the term corresponding to £ = 3™ is 1/(3™ + 1) and it 
does not contribute to 6*1. □ 

Step 2. If i(3™ - 1) < fc < 3™ - 1, then Si{n, k) is the empty sum. Therefore the 
denominator of hn,k is not divisible by 3™. 
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Proof. The index £ in the sum defining hn,k satisfies 1 < £ < 2 ■ 3"' ~ 1 — 2k. The 
assumption on k guarantees that neither ^ = 3™ nor £ = 3"^ — 1 appear in this 
range. □ 

Step 3. If fe = ^(3™ — 1), then the denominator of Si{n, k) is not divisible by 3"*. 

Proof. In this case the sum S\ (n, k) is 

1 _^3'" + l_3'" + 3 ^ 



3™(3™ + 1) (3™ - 1)3™ 32' 

Step 4. If < fc < 5(3™ — 1), then the denominator of Si{n, k) is not divisible by 
3". 

Proof. The proof of this step employs a theorem of Kummer stating that fp((^)) 
is equal to the number of borrows involved in subtracting h from a in base p. By 
Kummer's theorem, {2k+i) (2fe+i) divisible by 3, so neither of the two 
terms in S\{n, k) has denominator divisible by 3™. □ 

Step 5. If A; = 0, then the denominator of 5i(n, k) is not divisible by 3™. 

Proof. For A; = we have 

and the two terms in Hn whose denominators are divisible by 3™ add up to 
(5.23) + 2 _ gjjj 2 . 3TO— 1 

with denominator not divisible by 3"*. □ 

It follows that, for n = 2 • 3™, the denominator of the term hn,k is not divisible 

by 3". Thus, i/3(7n) < m - 1. 

Claim: ^^{"in) > m — 1. This is established by checking that 3™~^ divides the de- 
nominator of hnfi- Indeed, there are six terms in hnfi = n—Hn whose denominators 
are divisible by 3™"^, and their sum is 

(5 24) y-^_ = _g^= 49 

1=1 

Therefore 3™^^ divides the denominator of /i„,o. This completes Case 1. 

Case 2. Assume now that n is not of the form 2-3'". This states that the base 3 
representation of n is not of the form 200 • • • OO3. 

Let r = [logg n\ , so that y is the largest power of 3 less than or equal to n. We 
show that f3(7n) = r by exhibiting a value of the index k so that the denominator 
of hn,k is divisible by 3*^. 

Step 1. Assume first that the base 3 representation of n begins with 1. Then 
choose /c = 0. As before, hnfi = n — Hn. Observe that each term in the sum 

(5.25) lcmil,2,...,n).Hn = ±'^^^^^^ 

is an integer. The condition on the base 3 representation of n guarantees that only 
one of these integers, namely the one corresponding to ^ = 3'', is not divisible by 
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3. Thus there is no extra cancehation of powers of 3 in and as a result the 
denominator of iJ„ is divisible by 3''. 

Step 2. Assume now that the base 3 representation of n begins with 2. Choose 
k = i(3'' + 3"^^"^). As in the discussion in Case 1, there are at most two terms in 
the sum 

" -l- 2fc / _ 

with denominator divisible by 3''. The sum of these terms is 

1 /n-3''\ 1 /n-3'' + l 



3'-(3'' + l) V2fc + iy ' (3''-l)3'"V 2fc + l 
If n = 2 mod 3, then Kummer's theorem shows that 3 divides the second binomial 
coefficient but not the first; otherwise, 3 divides the first binomial coefficient but 
not the second. Therefore hn,k has precisely one term with denominator divisible 
by 3*^. The argument is complete. □ 

Corollary 5.10. The 3-adic valuation of the denominator an,2 oi An,2{x) is 

^ --^ 1 ^sC*^' lcm(l, 2, . . . , n)) — 1 if n = 2 • 3™ for some m> 1 
1 1/3(71,! lcm(l, 2, . . . , ri)) otherwise. 

Note 5.11. It easily follows that the 2-adic valuation of 7„ is 

(5.28) i/2(7n) = i^2(lcm(l, 2, . . . , n)) - 1. 

Indeed, from (|5.18p one has J/2(7ri) < z/2(lcm(l, 2, . . . , n)) — 1. On the other hand, 
the denominator of /i„^o = n — Hn is divisible by the highest power of 2, i.e., by 
2Liog2"J^ which implies 

The proof of Conjecture 15.51 has been reduced to the identity 

(5.29) z/p(7„) = i/p(lcm(l,2,...,n)) 
for all primes p > 3. 

The sequence 2 • 3™ appearing in the previous discussion also appears in relation 
with the denominators of the harmonic numbers _ff„. As before, write 

(5.30) Hr, = JV- 

in reduced form. The next result considers a special case of the quotient Dn^i/Dn of 
denominators of consecutive harmonic numbers. The general case will be described 
elsewhere [23^. 

Theorem 5.12. Let n e N. Then D2-3--1 = 3^2.3-.. 

Proof. An elementary argument shows that i'2{Dn) = \^og2n\. Therefore iV„ is 
odd and I?„ is even. 
Observe that 

iV2.3..-l , 1 



(5.31) 



D2.3" £'2.3"-l 2 • 3" 

^ 2 • 3"iV2. 311-1 + Z32-3"-i 

" 2-3"i?2.3"-l 
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Therefore the denommator £'2-3" is obtained from 2 • 3"£'2-3"-i by cancehng the 
factor 

(5.32) w = gcd (2 • 3" • N2.3^^i + i?2.3'.-i, 2 • 3" • i?2.3"-i) • 
That is, 

(5.33) 2-3"-D2.3"-i =u.•£'2.3"• 
Lemma 5.13. The number w has the form 2" • 3*^, for some a,/3 >Q. 
Proof. Any prime factor p of w divides 

2 • 3" • (2 • 3" • N2.3n_, + i?2.3"-i) - 2 • 3" • Z?2.3"-i = 2^ • 3^" • 

Then p is a common divisor of 2 • 3" • iV2.3"-i and 2 • 3" • D2-3"-i. The harmonic 
numbers are in reduced form, so p must be 2 or 3. □ 

The relation (|5.33p becomes 2 • 3"Z32-3"-i = 2" • 3'^-D2-3", and replacing this 
in ([OT|) yields 

(5.34) 2" • 3^7V2.3" = 2 • 3"7V2.3'^-i + ^2.3"-i. 

Define t = [log2(2 • 3" - 1)J > 1 and write L'2.3'«-i = 2*6*2. 3"-i with C2.3'.-i an 
odd integer. Then (|5.34p becomes 

(5.35) 2"-i • 3^7V2.3'. - 2*-iC2.3"-i = 3"7V2.3"-i. 

A simple analysis of the parity of each term in (|5.35p shows that the only possibility 
is a = 1. 

The relation (|5.33p now becomes 

(5.36) 3"-i?2.3"-i =3''-i?2.3". 

In the computation of the denominator -D2-3" we have the sum 

(5.37) 1 + 1 + .^ 



2 3 3" 2 • 3" - 1 2 • 3" 

so that the maximum power of 3 that appears in a denominator forming the 
sum ([QT)) is 3". Simply observe that 3"+^ > 2 • 3" - 1. The combination of 
all the fractions in the sum (|5.37p with denominator 3" is 

(5.38) ' ' 2 + 1 ' 



3" 2-3" 2-3" 2-3"-i' 
It follows that the maximum power of 3 in (|5.37l) is at most 3"~^. 

The terms in (|5.37p that contain exactly 3"^^ in the denominator are 

(^■39) 3^' 2 • 3""!' 4 •3"-!' 5- 3"-!' 

and these combine with the two terms with denominator exactly divisible by 3" to 
produce 

/ 1 1 1\ 1 1 49 

(1+2 + 4 + 5; -3^+2:3^= 2073^- 

The rest of the terms in (|5.37p have at most a power of 3"^^ in the denominator. 
The total sum can be written as 

.g^^s 49 Xn ^ 49yn + 60Xn 
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and no cancellation occurs. Therefore 3"~^ is the 3-adic valuation of Write 
D2-3" — 3"^^ • i?2-3", where -E2-3'« is not divisible by 3. 
Now consider the denominator Z32-3"-i- Observe that 

^ ' ^ 2 + ■ ■ ■ + 2 • 3" - 1 " 3« ^ y„ • 3"-i " yn • 3" ' 

with yn not divisible by 3. Therefore 3" is the 3-adic valuation of -D2-3"-i- Write 
-D2-3"-! = 3" • £'2-3"-i where iJ2-3"-i is not divisible by 3. 

The relation ([06| now reads i"^"- E2.3^^-l = 3'^+"-i£'2.3" and this gives /3 = n+1. 
Replacing in (|5.36p produces D2-3"-i — 3L'2-3", as claimed. □ 

6. The iterated integral of ln(l + x^) 
In this final section we consider the iterated integral of ln(l + a;'^). The first value 



6 

- \/3arctan \^ ^ j + (x + I) \n{x + 1) + -{2x - l)In(a::^ - x+1). 
This and additional values suggest the ansatz 

(6.1) fn{x) = An,3{x) + Bn,3{x)u + Cn,3ix)v + Dn^^{x)w, 

where 

^ (\-2x^ 
M = V 3 arctan 



V3 

V — ln(x + 1) 

w — ln(x^ — x + 1) 

and where A„^3 is a polynomial in Q[A/37r, a;] and -Bn.s, C'„,3, and I?n,3 are polyno- 
mials in Q[x]. 

The method of roots described in Section [3] shows that /n(x) can be expressed 
in terms of ln(a; + 1), ln(a; -I- w), and ln(a; -I- w), where cj = g2i-»/3 = ^(—1 -I- iv^) 
satisfies ufi — \. Using the relation 

1 / 1 — 2.x \ 

(6.2) ln(-2i(a; + w)) = - ln(x^ - a; + 1) + i arctan \^ ^ j + ln(2), 

we can convert between (j6.ip and expressions in terms of ln(a:; + w) and ln(a; -I- w). 

As was the case for the iterated integrals of ln(l -I- x) and ln(l + a;^), it is easy 
to conjecture closed forms for all but one of these polynomials. 



Theorem 6.1. Define 



if fc = mod 3 
X3(fc) = 1 iffc = l mod 3 
-1 if fc = 2 mod 3 



and 

A(fc) 



1 if fc = mod 3 
if fc ^ mod 3. 
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Then 

1 " / 



fc=0 ^ 



fc=0 ^ 

Proof. The method of roots developed in Section [3] shows that the iterated integral 
can be expressed in the form (|6.ip . The polynomials An^s, Bn,3,Cn.3, Dn.3 will be 
linear combinations of the powers (a; + 1)", (x + w)", and (a; + w)"- Comparing 
initial values, it is found that 

i?„,3 (x) = -^((x (a: +(^r) 
c„.3(x) = -Ax+ir 

n] 

^„,3(x) = ^((x+c.r + (x+c^)"). 

Note that the above expressions can also be automatically found as solutions of the 
fourth-order recurrence that HolonomicFunctions derives in this case: 

(6.3) {n-2){n~l)n^Fn 

= (n - 2)(n - l)(4n - 3)xFn-i - 3{n - 2)(2n - 3)x'^ Fn-2 

+ [{4n - 9)x^ + n] F„_3 - x{x^ + 1)F„_4. 

The above closed forms for -6,1,3 and £'n,3 can be used to derive explicit expres- 
sions for their coefhcicnts: 

The value of the last parenthesis can be found by case distinction using the fact 
that oj^ = u)^ = 1: 



n — k = mod 3 
n — k = I mod 3 
n — k = 2 mod 3 



1-1 = 0, 

U! — U! — iVS, 
;2 _ ^2 ^ 



It follows that 



' fc=0 ^ ^ 



x\ 



The similar computation for £',1,3(0;) is left to the reader. □ 

Note 6.2. In order to obtain these functions from a purely symbolic approach, 
consider a brute force evaluation of fn{x) by using Mathematica to evaluate (|1.4p . 
The results are expressed in terms of the functions 

(6.4) h,{x) = 2F,[ ^(^'3^ and /i2(x)=2Fi( ^^'3^ ; -rr^) , 
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where 

(6.5) ^^i(";';-)=E 



a, 6 _ (a)fc {b)k^k 



(c)fcfc! 

is the classical hypergeometric series. The first few values are 
fi{x) = -3x + X ln(l + x^) + 3xhi{x) 

/2(x) = + ln(l + x^) + 3x^/11 (x) - ^h2{x) 

Mx) = + lix^ + l)Hl + x') + ^-^h,{x) - ^-^h,{x). 

The hypergeometric series hi {x) and /12 (a;) can be expressed as 

hi {x) — — (ln(l + a;) + w ln(l + ujx) + w ln(l + ux)) 
ix 

2 

^12(2;) ~ 7—^ (— ln(l + x) — uj ln(l + cjx) — UJ ln(l + wx)) , 

with w = e^'^*/^ as before. These expressions can be transformed into the functions 
in (ICT . 

As in previous sections, the closed-form expression for the pure polynomial part 
An,z{x) is more elaborate. By writing 

fn{x) = An,z{x) + Bn.six) \n{x + 1) + Cn^six) \n{x +Uj) + DnA^) + 

we obtain a polynomial An ^{x) with rational coefficients. The first values are given 

by 

(6.6) Ao,3{x) = 0, ii,3(x) = -3x, A2,3{x) = -fx^, A^.^ix) = -^x\ 

Schneider's Mathematica package Sigma can be used again to obtain, from the 
recurrence (|6.3p and the initial conditions, the expression 

1 " fe 

(6.7) A„,,ix) ^ V ^ [(a: + l)""'' + {x + 0.)"-'= + (x + . 

k=l 

6.1. Arithmetical properties of An^s- Define a„_3 to be the denominator of v4„_3 
and /3„,3 = a„,3/(na„_i,3). 

Conjecture 6.3. The sequence /3„_3 is given by 

p if n = ^ 3 for some prime p and m G N 
if = 3 • 11'" for some m e N 



(6.8) 



11 

11 if n = 3 • 11'" + 1 for some m e N 
1 otherwise. 



Observe that this expression for /3„_3 does not have the exceptional case where 
3 ■ 11™ + 1 is a prime power that appears in /3„^2 given in (|5.8p . This is ruled out 
by the following. 

Lemma 6.4. Let m G N. Then 3 • ll'" + 1 is not a prime power. 

Proof. The number 3 • ll'" + 1 is even, so only the prime 2 needs to be checked. 
We have 3 • 11™ + 1 = 3'"+^ + 1^0 mod 8 since the powers of 3 are 1 or 3 modulo 
8. Therefore 3 • 11™ + 1 (since it is larger than 4) is not a power of 2. □ 
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